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We present a split-beam neutron interfero- 
metric experiment to test the non-cyclic 
geometric phase tied to the spatial 
evolution of the system: the subjacent 
two-dimensional Hilbert space is spanned 
by the two possible paths in the interfero- 
meter and the evolution of the state is 
controlled by phase shifters and absorbers. 



A related experiment was reported 
previously by some of the authors to 
verify the cyclic spatial geometric phase. 
The interpretation of this experiment, 
namely to ascribe a geometric phase to 
this particular state evolution, has met 
severe criticism. The extension to non- 
cyclic evolution manifests the correctness 
of the interpretation of the previous 
experiment by means of an exphcit 
calculation of the non-cychc geometric 
phase m terms of paths on the Bloch- 
sphere. The theoretical treatment 
comprises the cyclic geometric phase 
as a special case, which is confirmed by 
experiment. 
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1. Introduction 

Since the discovery of a geometric effect by Berry 
[1] in the shape of an additional phase factor after an 
adiabatic and cyclic transport of a quantum system, 
Berry's phase has been intensively investigated and 
generalized: the extension to degenerate subspaces by 
Wilckzek [2], the removal of the adiabatic constraint by 
Aharonov and Anandan [3] and the cyclic condition by 
Samuel and Bhandari [4] using the early ideas of 
Pancharatnam [5] and the kinematic approach to 
geometric phases by Mukunda and Simon [6], to name 
a few. In all theses contexts the geometric phase is 
dependent only on the geometry of the subjacent 



Hilbert space, but not on the particular dynamics of the 
system under consideration. Furthermore, Manini and 
Pistolesi [7] proposed an off-diagonal geometric phase 
to exhibit the geometry of state space in situations 
where the usual (diagonal) geometric phase is unde- 
fined. This has been verified experimentally by some of 
the authors [8] . 

In the course of the development of quantum 
mechanics it has become clear that the concept of pure 
states is not sufficient when taking environmental influ- 
ences causing decoherence effects into account. Then 
one has to use the concept of mixed states. Probably the 
first treatise of a geometric phase for mixed states is 
due to Uhlmann [9] in a quantum algebraic context. 
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Another definition of a mixed state geometric phase is 
due to Sjoqvist et al. [10] using an interferometric 
approach for its definition. For these concepts one has 
to keep in mind that there exist points in parameter 
space for which the mixed state geometric phases 
remain undefined provoking an extension to off- 
diagonal mixed state geometric phases [11]. 

The geometric phase is associated with an evolu- 
tion of a system governed by an Hamiltonian, e. g., a 
neutron in a magnetic field where the geometric phase 
arises by the spinor evolution due to the coupling with 
the magnetic field. Here, we observe a geometric phase 
as an effect of the change in the spatial degrees of free- 
dom in an interferometry setup. A proposal to verify the 
spatial geometric phase is due to Sjoqvist [12] using 
polarized neutrons by reversing the roles of the magnet- 
ic field and the spatial degrees of freedom. Moreover, 
an experiment using unpolarized neutrons has been 
performed by Hasegawa et al. [13, 14] to test the cychc 
spatial geometric phase by inducing a relative phase 
shift of 2n between the interfering neutron beams in a 
perfect silicon single-crystal interferometer. The geo- 
metric interpretation of this experiment has been 
dismissed by Wagh [15] demanding further investiga- 
tions, namely in the non-cyclic case, which is the 
purpose of the current article. 

2. Geometric Phases 

Let us briefly review the basic concepts of geometric 
phases: A geometric phase is a quantity which is 
deeply connected to the curvature of some underlying 
(state- or parameter-) space. A two-dimensional plane 
in three-dimensional real space does not have an 
intrinsic curvature, but when considering a sphere 
embedded in euclidean real space, we have to take the 
curvature of this manifold into account. In geometry 
this curvature is reflected, for example, in the angle 
difference of a vector transported around a loop along 
geodesies, i. e., great circles: If a vector is pinned onto 
a sphere and then transported along a meridian to the 
equator, for some angle a along the equator and back to 
the initial point without changing its length and its 
direction in the tangent plane to the surface of the 
sphere, the vector will point in a different direction with 
a relative angle of a as the holonomy associated with 
the loop. If we do the same on a two-dimensional plane 
the initial and the final vector will point in the same 
direction. 

Berry [1] was the first who addressed this issue in 
quantum mechanics: He considered a system initi- 
ally in an eigenstate \n(R(t)))i=Q of the governing 



Hamiltonian H{R{t)) dependent on the parameters R{t) 
changing with time /. As a demonstrative example one 
may consider a neutron coupling to a magnetic field 
H(R(t)) = -ju ' B(R(t)) due to its magnetic dipole 
moment ju = ju^G, where 0"= {cr^, o^, o;} are the Pauli 
matrices and ju^ denotes the magnetic moment of a neu- 
tron. Suppose now that the neutron is initially polarized 
in the direction of the magnetic field. If the direction of 
the magnetic field is changed adiabatically, i. e., slow- 
ly enough to avoid transitions to an orthogonal state, 
the system will stay in the eigenstate \n(R(t))} at all 
times /. Furthermore, when tracing out a loop in para- 
meter space the final state | *f^(T)) at time T will be the 
same as the initial state up to an additional phase factor: 



|»F(T)) = e'^V'^^|n(/f(T))) 

— i^E^iOdt -S {n(R)\V^\n(R))dR 



n(R(T))). 



(1) 



1 rr 

The first phase value d= ~\ K (^(O)d^ is dependent 

on the time needed to traverse the loop and on the in- 
stantaneous energy E^(t) = {W(t)\H(t)\ W(t)} of the system, 

whereas the second phase (l>g = i% (ri(R)\V j^\n(R)) 

dR is dependent only on the circuit integral in para- 
meter space revealing the geometric structure. The latter 
is termed Berry phase or more general geometric phase 
in contrast to the former dynamical phase i^j. <^g can 
be rewritten as a surface integral by use of Stoke 's 

Theorem yielding <^^ =-lmj dSV^{RX where F is 

the surface enclosed by the loop in parameter space 
with dS denoting the area element and V^= Vx{n\ Vn) 
in an obvious abbreviated notation. For the neutron 
example — or more generally for any spin- 1/2 
particle — (j)^ equals half of the solid angle enclosed by 
the loop as seen from the degeneracy point |i?| = in 
parameter space. This can also be related to the 
example from geometry above where the holonomy 
after the transport of the vector pinned initially to the 
north pole of a sphere equals the solid angle as seen 
from the origin of the sphere. 

Several restrictions have been relaxed in course of 
the years, e. g., extensions to nonadiabatic [3], noncyclic 
and nonunitary [4], and nonpure [9, 10] geometric 
phases have been made. Important in our case are the 
generalizations to the nonadiabatic regime by 
Aharonov and Anandan and to noncyclic paths by 
Samuel and Bhandari. 

In this case we have to introduce the Projective 
Hilbert space (Ray space) H by identifying all state 



252 



Volume 110, Number 3, May- June 2005 

Journal of Research of the National Institute of Standards and Technology 



vectors in Hilbert space J^ which differ only by an 
overall phase factor: 



|f) = e'«|<^),a( 



(2) 



The stress is therefore shifted from the parameter space 
of the Hamiltonian in case of Berry's construction to 
state space. We are not interested in the changes of the 
driving parameters (as the direction and strength of the 
magnetic field) but in the changes of the state itself In 
Berry's considerations these two spaces are identical 
since the state follows the changes in parameters due to 
the adiabaticity condition. 

In the construction by Aharonov and Anandan one 
considers an open path in Hilbert space which is 
projected to a path in Ray space by use of the equiva- 
lence relation in Eq. (2), i. e., the curve C : / e [0, r] ^ 
{(j) (t)) e H is projected to C: ^ e [0, r] -> 7i;(|<^ {t))) = 
\(l) (t)) {(j) it)\ e n with |<^(0))~|<^(T)). For Can 
absolute phase factor of |<^(0) is immaterial, since the 
curve C is defined via the evolution of the projection 
operator |<^ (t)) {(j) {t)\. The geometric phase is a proper- 
ty of Ray space, where C is closed due to the equiva- 
lence of the initial and the fmal state |<^(0))~ |<^(t)) 
and can be calculated via a surface integral over the 
area enclosed by C. 

One can fmd many different curves C, C\ ... in 
Hilbert space differing by a phase factor e'^^^ and yield- 
ing the same curve in Ray space under the projection 
map K. On the other hand for a given curve in Ray 
space, there exists one distinct curve in Hilbert space 
fulfilling the parallel transport conditions, namely that 
two neighbouring states |<^(0) and |<^(/ + dO) in H 
have the same phase, that is to say, (<^ (0|<^ (t + d^) is 
real and positive. This implies by Taylor expansion that 

(0(/)| — |0(/)) = O[18]. For this curve the dynamical 

phase vanishes as one can verify by inserting the 
Schrodinger equation 77 (^)|0(^)) = /M/d^|0(/')) into the 
parallel transport condition. 

The concept can be extended to apply to open paths 
in Ray space where |<^ (r)) -/- 1<^ (0)) by closing the curve 
by a geodesic, i. e., a path in Ray space with the short- 
est distance from |<^ (r)) <<^ (t)| to |<^(0)) <<^(0)|. Then 
one obtains a well-defined surface area enclosed by the 
path generated by the evolution of the system plus the 
geodesic closure. This surface provides an expression 
for the geometric phase, which has been proven by 
Samuel and Bhandari [4]. 

To sum up, for a general evolution of a quantum state 
the state obtains a dynamical phase dependent on the 



energy and time as well as a geometric phase only 
dependent on the subjacent geometry of state space. 
For special Hamiltonians which fulfill the parallel 
transport conditions the dynamical phase vanishes, 
which is also the case when the state is transported 
along a geodesic. An example of the latter is an evolu- 
tion along a great circle on a sphere for a two-level 
system which we will encounter in the forthcoming 
discussion. 

3. Interferometric Setup 

Due to Feynman [ 1 6] the description of any two- 
level quantum system is equivalent to the description of 
a spin- 1/2 particle. Exploiting this equivalence there is 
in principle no difference between manipulations in 
the spin space of neutrons with the orthogonal basis 
{|T)^ \i)} as eigenstates of c^, and momentum space 
with {| ^ ), I k' )} as orthogonal basis vectors correspon- 
ding to two directions of the neutron beam in an inter- 
ferometer. In both cases one can assign a geometric 
phase to the particular evolution of the initial state. An 
even more appropriate description for the interfero- 
metric case for the forthcoming discussion is in terms 
of "which- way" basis states {\p),\P^)}^ namely, if the 
neutron is found in the upper beam path after a beam- 
splitting plate it is said to be in the state \p), or in the 
state I p "^), if found in the lower beam path. In case of 
a 50 : 50 beamsplitting of the incident (neutron) beam 
into a transmitted beam and a reflected beam, the asso- 
ciated wave vector after the beamsplitter can be written 
as an equally weighted coherent superposition of 
the two paths | ^(5) ) = l/^(\p)-\-&^\p'-)^with the 

relative phase 5 e M depending on the particular phys- 
ical realization of the beamsplitter. 

For testing the spatial geometric phase we use a 
double-loop interferometer (Fig. 1), where the incident 
unpolarized neutron beam | y/) is split up into a diffracted 




Fig. 1. Experimental setup to test the spatial geometric phase in a 
neutron interferometer. 
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reference beam | (//^ef) and a transmitted beam | y/^}. The 
transmitted beam is subjected to further evolution in 
the second loop of the interferometer by use of beam- 
splitters (BSl and BS2), an absorber (A) with trans- 
mission coeficient Tand the phase shifter PS2 generat- 
ing a phase shift of e'^^ on the upper and e'^^ on the 
lower beam path, respectively, yielding the final state 
\y/{) = ^ \¥t) ^ ^ \p)' H^re \y/t) = \p) since before the 
beam splitter BSl the beam is clearly localized as seen 
from the second loop of the interferometer. The unitary 
matrix U= U(T, (f)^, 1^2) comprises all the manipulations 
in the second loop: 



4~2 

1 



41 



(e'*» + Vre'*^|/>^))^f/|v/.) = |r,)- 



(3) 



The geometric phase can then be extracted from the 
argument of the complex valued scalar product 
between the initial and the final state arg {yf\yf^ (when 
removing dynamical contributions as will be discussed 
later). This is where the reference beam comes into 
play: | y/^^}^ is not subjected to any further evolution, but 
is stationary apart from adding a phase factor e'^ by use 
of the phaseshifter PSl. \\ff,^^ propagates towards the 
beamsplitter BS2 from the upper path, thus we can 
assert it to be in the state e'^|/?). Then by the variable 
phase shift e'^ one can measure the shift of the interfer- 
ence fringes reflecting the phase difference between 
IV^Jandlv^^. 

This preparation of the states is followed by the 
recombination of the two beams \\ff}) and \\ff^^ at the 
beamsplitter BS2 and the detection at the detector D^ in 
the forward beam. This step can be described by the 
apphcation of the projection operator \q) {q\ = l/2(\p) + 
\P^))({p\'^{p'\) (with 5=0, which can always be 
achieved by an appropriate choice of the phase of the 
basis states) to ly/^ as well as to ly/mf)- 






(4) 



where K is some scaling constant. 

The intensity / measured in the detector Dq is pro- 
portional to the absolute square of the superposition 

\¥'f) + ^'WJ' 

+2|Kf |v^;)|cos(77-arg(i^,:f |i^;)). (5) 

We notice a phase shift of the interference pattern 
by arg (y/'mflY'f)- This phase shifts corresponds to the 



Pancharatnam connection [5] between the state 
ly/'f) and the state ly/'J) = |^) (^Iv^t) = k) from which 
we can extract the geometric phase. Explicitly we 
obtain 

<^ = arg«f|v/^;) 

01+02 



- arctan 



-if 



1-Vr 



i+Vr 



(6) 



where A<^= 02~0i- The geometric phase is defined as 
[6] 



<^g=arg(v^;;^|v^;)-^^, 



(7) 



where (j)^ denotes the dynamical part. From Refs. [13] 
and [15] we know that the dynamical part stemming 
from the phase shifter PS2 is given by a weighted sum 
of the phase shifts <^i and <^2 with the weights depend- 
ing on the transmission coeficient T. In particular we 
have 

^^=-i7^' (^) 

which vanishes by an appropriate choice of phase shifts 
and transmission, i. e., (1)^ = for i^^ /<^2 = -^^ 

By varying the relative phase A<^ from to 2n and 
setting 0d = the geometric phase <^g can be plotted 
over A<^ (Fig. 2). 
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Fig. 2. Geometric phase ^„ [rad] in dependence on the relative phase 
shift A(^[rad] for r= 1/8. 



4. Bloch-Sphere Description 

For every two-level system we can use the Block- 
sphere for depicting the state vectors and evolutions 
thereof as points and curves on a sphere. Then the 
results obtained above, i. e., the shift of the interference 
pattern in Eq. (5) without dynamical contributions, 
should be equal to the (oriented) surface area enclosed 
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by the paths of the state vectors on the Bloch-sphere, or, 
equivalently to the solid angle traced out by the state 
vectors as seem from the origin of the sphere. 

As we can observe in Fig. 3, the north pole of the 
sphere can be identified with a state with well known 
path, i. e., an eigenstate of the observable |/?) {p\. After 
the beam splitter BSl the state |i//t) evolves to an equal 
superposition of upper path and lower path, there- 
fore the evolution on the Bloch sphere is given by a 
geodesic from the north pole to the equatorial line (the 
particular point on the equator is arbitrary due to the 
arbitrary choice of the phases of the basis vectors). 




Fig. 3. Path of the state in an interferometer on the Bloch sphere 
representing the 2-level system (upper path \p) {p\ and lower path 

\p')(p'\y 



The absorber changes the weights of the superposed 
basis states, in particular for the extremal values of T 
parameterized by the angle with r= tan^ 0/2, we end 
up either again with an equally weighted superposition 
for no absorption (T=l or = nil) or the state is now 
on the north pole for total absorption (7= or ^ = 0), 
since in the latter scenario we know the particle has 



taken the upper path when detecting a neutron in Dq. 
For Te (0, 1) the state is encoded as a point on the 
geodesic from the north pole to the equatorial line. 

Due to the phase shifter PS2 we obtain a relative phase 
shift between the superposing states of A<^= (l>2~^\' 



;^(I^)H^^»-;^(I^)+^'"I^^» (9) 



This can be depicted as an evolution along a circle of 
latitude on the Bloch sphere with periodicity of 2;:. 

The recombination at BS2 followed by the detection 
of the forward beam in Dq is represented as a projection 
to the starting point on the equatorial line, i. e., we have 
to close the curve associated with the evolution of the 
state by a geodesic to the point \q) {q\ on the sphere as 
discussed for non-cychc paths in Sec. 2. As for the 
reference state | i/^ref) we note that the phase shift of t] 
has no impact on the position of the state on the Bloch 
sphere, it stays at the north pole. Due to the recombina- 
tion at BS2 and the detection the state is also projected 
to 1^) (^1 contributing to the forward beam incident to 
the detector Dq. 

The paths are depicted in Fig. 4 in detail for cyclic 
Fig. 4(a) as well as non-cyclic evolution. Fig. 4(b). For 
a relative phase difference greater than Till we have to 
take the direction of the loops into account. In Fig. 4(b) 
the first loop is traversed clockwise, whereas the 
second loop is traversed counterclockwise yielding a 
positive or negative contribution to the geometric 
phase, respectively. 

The numerical calculation of the surface area F 
enclosed by the path traversed by the neutron is straight- 
forward by evaluating the solid angle Q. = \f sin d(A<^)d^ 
via a surface integral and using <^g = -Q./2. For the cyclic 
case this integral can be solved easily by calculating the 
segment on the sphere according to Fig. 4(a) to obtain 
<^g = -Q/2 = 7r(cos^-l). For the non-cyclic case no 
analytic expression has been found to compare the 
results with the phase shift of the interference fringes 
appearing in Eq. (5). However, the numerical results 
are equivalent and agree with Fig. 2. 
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(a) Cyclic evolution (b) Non-cyclic evolution 

Fig. 4. Paths on the Bloch sphere corresponding to the evolution of the state in the splitbeam experiment. 



5. Experimental Results 

For an experimental test of the spatial geometric 
beam we have used the double-loop perfect-crystal- 
interferometer installed at the S18-beamline at the 
high-flux reactor at the Institut Laue-Langevin, 
Grenoble [ 1 7] . A schematic view of the setup is shown 
in Fig. 1. Before falling onto the skew-symmetric 
interferometer the incident neutron beam was collimat- 
ed and monochromatized by the 220-Bragg reflection 
of a Si perfect crystal monochromator placed in the 
thermal neutron guide H25. The wavelength was tuned 
to give a mean value of /Lq = 0.2715 nm. To eliminate 
the higher harmonics we have used prism-shaped 
silicon wedges. The beam cross-section was confined 
to (5 X 5) mm^ and an isothermal box enclosed the 
interferometer to achieve reasonable thermal environ- 
mental isolation. For the phase shifters parallel sided 
aluminium plates have been used: 4 mm inserted in the 
first loop as PSl and 4 mm and 0.5 mm, respectively in 
the second loop for PS2 yielding a ratio of 1/8 for <^i/<^2- 
To avoid dynamical phase contributions a gadolinium 
solution with 7=0.118 [19] has been used as an ab- 
sorber. For a comparison with the theoretically predict- 
ed values one has to keep in mind that the contrast 
reflecting the coherence properties is different between 
each of the beams in the interferometer. Accounting for 
this experimental fact in the theoretical derivation of 
the geometric phase we notice a slightly flattened 
curve in Fig. 5 compared to Fig. 2. Nevertheless, one 
can recognize the increase in geometric phase for 
Ai^e [0, n/2] due to the positively oriented surface 
followed by a decrease due to appearance of a counter- 



clockwise traversed loop on the sphere yielding a 
negative phase contribution. 
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Fig. 5. Experimental verification of the spatial geometric phase 
using a neutron interferometry setup. 



6. Conclusions 

In summary we have shown that one can ascribe a 
geometric phase not only to spin evolutions of neu- 
trons, but also to evolutions in the spatial degrees of 
freedom of neutrons in an interferometric setup. This 
equivalence is evident from the description of both 
cases via state vectors in a two dimensional Hilbert 
space. However, there have been arguments contra to 
the experimental verification in [13] which we believe 
can be settled in favour of a geometric phase appearing 
in the setup described above. The twofold calculations 
of the geometric either in terms of a shift in the inter- 
ference fringes or via surface integrals in an abstract 
state space allows for a geometric interpretation of the 
obtained phase shift. 
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